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ABSTRACT 

Current measurements of the large-scale structure of galaxies are able to place an 
~ 0.5 eV upper limit on the absolute mass scale of neutrinos. An ordcr-of-magnitude 
improvement in raw sensitivity, together with an insensitivity to systematic effects, is 
needed to reach the lowest value allowed by particle physics experiments. We consider 
the prospects of determining both the neutrino mass scale and the number of of massive 
neutrinos using future redshift surveys, specifically those undertaken with the Square 
Kilometre Array (SKA), with and without additional constraints from the upcoming 
Planck CMB experiment. If the sum of the neutrino masses X] > 0.25 cV then the 
imprint of neutrinos on large-scale structure (LSS) should be enough, on its own, to 
establish the neutrino mass scale and, considered alongside CMB constraints, it will 
also determine the number of massive neutrinos N^^ and hence the mass hierarchy. 
If ^ mi 0.05 eV, at the bottom end of the allowed range, then a combination of 
LSS, CMB and particle physics constraints should be able to determine '}^mi, N„ 
and the hierarchy. If is in the specific range 0.1 — 0.25 eV, then a combination 

of LSS, CMB and particle physics experiments should determine "^mi, but not 
or the hierarchy. Once an SKA-like LSS survey is available there are good prospects 
of obtaining a full understanding of the conventional neutrino sector, and a chance of 
finding evidence for sterile neutrinos. 

Key words: Cosmology: cosmological parameters - Cosmology: cosmic microwave 
background - Cosmology: large-scale structure of the universe - surveys - neutrinos 



1 INTRODUCTION 



The proof that neutrinos have mass was a major break- 
through in particle physics. Th is proof came from the obser- 
vation (e.g. lFukuda et al]ll996l ) that neutrinos in one weak- 
flavour state are able to 'oscillate into' neutrinos of a differ- 
ent weak-flavour state which, in 'the vacuum', is disallowed 
by quantum mechanics unless neutrinos have mass. The im- 
plications of this result are profound not only for particle 
physics, but also for cosmology. It was the first measured 
effect that is not included in the Standard Model of parti- 
cle physics. It also started to set constraints on the abso- 
lute meiss scale of a particle species which cosmologists have 
yet to include properly in the energy density budget of the 
Universe. A r ecent review of the critical issues is given by 
iKavseil (|2005l ') and another recent review approaching neu- 
trino physics from the theore tical cosmology side is given in 
iLesgourgues fc Pastod (|2006l '). 

Throughout this paper we keep with particle physics 
notation by using natural units in which c = h = k — 1. 
We will refer throughout to a 'fiducial' cosmology in which 



the normal cosmological constants take the values: {Qt, fic, 
h, Us, erg} = {0.04, 0.26, 0.72, 1.0, 0.9}. We define Q.m as 
the fraction of critical density contributed by all matter: 
baryons, CDM and neutrinos. The power spect rum P(k) is 
calcu lated using the 'Boltzmann code' CAMB (|Lewis et al.l 
l2000h . We reserve the use of the symbol rui, for the absolute 
mass scale of neutrinos, by which we mean the rest mass 
of the most massive neutrino (see Fig. [T|. Throughout this 
paper we assume that the Universe is flat and that the 'dark 
energy' is Einstein's cosmological constant with an equation- 
of-state parameter w = 



1.1 Background Particle Physics 

Vacuum neutrino oscillations occur if neutrinos are massive 
because their mass eigenstates are an admixture of the weak- 
flavour eigenstates familiar from the physics of electroweak 
interactions. This can be written as 
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Figure 1. The (mass)'^ spectrum of neutrinos allowed by current 
neutrino oscillation experiments for (left) a 'normal' hierarchy 
and (right) an 'inverted' hierarchy. We define the absolute mass 
scale of the neutrino rrii, as the mass of the most massive cigcn- 
state. This scale is bounded below by the minimum required by 
oscillation experiments (~ 0.05 eV) and above (at ~ 0.5 eV) by 
current cosmological ( l arge-scale structure) measurements 
lElgarOv fc Lahavlliooa: ISeliak et allbOOSl: iFukugita et al.ll200l 
some of which are expected to have more systematic effects than 
others. If the true value of lies towards the top of this allowed 
range, then we call this a quasi-degenerate scenario because, re- 
gardless of the hierarchy, mi m2 ~ ms (with m,^ = ma for 
a normal hierarchy, and = m2 for an inverted hierarchy). If 
the true value of rrti, is towards the bottom of the allowed range 
= 1713 (with mi m2 ~ 0) for a normal hierarchy; and 
= 1712 ~ mi (with ms 0) for an inverted hierarchy. The 
size of each coloured bar is proportional to the contribution of 
weak- flavour eigcnstates to each mass eigenstate: red, fe[|f/ciP]; 
green, u^l\U„i\'^]] and blue, i/rHC/riP]- 



(1) 



where a = 1,2,3 labels the weak-flavour eigenstates and 
i = 1, 2, 3 labels the mass eigenstates. Ud is the mixing 
matrix that governs the rate of oscillations between different 
weak-flavour eigenstates. This yields a probability of a given 
neutrino oscillating from a weak-flavour state a to another 
weak-flavour state /? of 

p.p = I {'^pWAt)) \" = I EE i'^MlMt)) f- (2) 



between muon and tau anti-neutrinos, controlled largely by 
Am32; experiments that probe electron anti-neutrinos orig- 
inating in the sun find strong evidence for a different, but 
again approximately two-neutrino, oscillation between elec- 
tron and muon neutrinos, controlled largely by Amf2. The 
current sta tus of quantita t ive ob servational results are sum- 
marised bv lMahoni et al.l (|2004l ) (see also Fig. [l]). 

The mixing matrix can be written in terms of a com- 
bination of Euler rotations in an ti = 3 dimensional space, 
namely 
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ing angles' 9ij. We have assumed that neutrinos behave like 
Dirac particles - so that each neutrino and its associated 
anti-neutrino are distinguishable - and we have introduced 
a phase 5 which allows for Charge-Parity (CP) violations 
in the sense of small differences between the probability of 
oscillations between neutrinos and 'CP-mirror-image' oscil- 
lations between anti-neutrinos. 

The transition probability can now be written as 



1 

4E 



Pc0 = S^p - 4:^^Re[UciU^iU*jU,3j] 

(4) 

where L is the distance between the source and the observer, 
which is assumed to be a vacuum, and Anvfj = rnl — 
because vacuum oscillations conserve energy, E = Ea = Ejs. 
From Eqn.|4]it is clear why particle physics experiments are 
mainly sensitive to the differences between the squares of 
the masses of the mass eigenstates and the mixing angles. 

The three main goals of future research into the neu- 
trino sector are as follows; (i) to determine the absolute mass 
scale for neutrinos; (ii) to determine whether the hierarchy is 
normal or inverted (see Fig. [1} ; (iii) to check whether there 
are so-called sterile neutrinos, particles which interact only 
via gravity, additional to the three weak-fiavour eigenstates. 
The prospects of rapid progress via fu rther particle physics 
experiments is limited (|Kavseij |2005| ) and strongly depen- 
dent on unknown factors like the true magnitude of 6'i3 , and 
whether or not the extra diagnostic power available, in prin- 
ciple, from measuring oscillations in the presence of matter 
can be successfully harnessed. 



The relationships between the weak-flavour and the 
mass eigenvectors for neutrinos are illustrated in Fig. [T] The 
third mass eigenstate is almost perpendicular to the axis 
of the electron neutrino eigenstate, or to put it in another 
way the 'mixing angle' defined by sin^ 613 — |(7c3p is small. 



The 



exp 



erimcntal evidence for oscillations is summarised 



by iKavser (,2005. ), but can be very crudely summarised as 
follows; experiments that probe muon and electron anti- 
neutrinos originating from cosmic ray showers in the atmo- 
sphere can be approximated by a two-neutrino oscillation 



1.2 Background Cosmology 

Neutrinos decoupled from radiation and matter within the 
first few minutes of the history of the Universe. As this 
de-coupling occurred before electron-positron annihilation 
produced an extra source of entropy for the photons, the 
neutrino temperature has remained at a fraction of the 
temperature T-^ of the CMB given by 
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(5) 



Neutrinos remain relativistic until the point at which 
~ m^, which corresponds to a redshift z ~ 6000 for 
rrii, ^ 1 eV or « ~ 30 for rrii, ~ 0.05 eV. This means that 
neutrino free streaming can suppress structure formation on 
small scales at early times, but the neutrino energy density 
is dominated by the rest mass at recent times. We assume 
that there are N^, 'massive' neutrinos where Ni, ~ 3 for 
a quasi-degenerate scenario, Ni, ~ 2 for an inverted, non- 
quasi-degenerate scenario and A*'^ ~ 1 for a normal, non- 
quasi-degenerate scenario (see Fig. [!}. We further assume 
that neutrinos are Dirac particles so that each neutrino and 
its associated anti-neutrino are distinguishable. We can then 
relate the energy density of photons (bosons) to the energy 
density of neutrinos (fermions) . The number density of neu- 
trinos can be related to the number density of photons by 
rii, = (3/ll)n-yA^^, so we can write 



3 m^N^n-, 



93.14/i2 eV 11 pc 



93.14/i2 eV 



(6) 



Neutrinos behave as a component that transitions from 
having an equation-of-state p = w^p, where p and p are the 
pressure and energy density respectively, with Wi, = 1/3 to 
— > as they become non-relativistic. We can evaluate 
(z) for neutrinos using 



il =0.05 N =1 



k h Mpc" 

Figure 2. The power spectrum P{k) for our fiducial cosmology 
(solid black line) compared with those with difi^erent values of 
and Af^, and through Eqn|6] (coloured dashed lines). There 
is a clear damping length scale imprinted on P{k) due to the 
free-streaming of neutrinos at early times (Eqn. |8j. This scale 
is obviously dependent on (Eqn. [9]l and more subtly on Ni, 
because, for a given C^, a smaller A^^ delivers a larger which 
goes non-relativistic earlier and hence produces less damping. LSS 
data should hence be able to constrain both the mass and the 
number of massive neutrino eigenstates. 
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where g is the degeneracy and / is the phase-space distribu- 
tion function for fermions. 

The effect of neutrinos on large-scale structure is illus- 
trated in Fig. [21 The main feature is a strong attenuation 
of the power spectrum P{k) on small scales if there is a sig- 
nificant energy density of neutrinos. This is associated with 
a preferred length (and k) scale given by the horizon when 
the neutrinos b ecome non-relativistic which, from Eqn. 1 of 
||Hu et aLlligoi ). is 



K ==; 0.026 ( ^^"l ^''^ ^11'^ h Mpc 
V 1 eV/ 



(8) 



assuming a quasi-degenerate scenario, where is the asso- 
ciated co-moving wavenumber. The damping of t he power 
spectrum at large scales, small k, is (Eqn. 2 of IHu et al.l 
1 19981 ) given approximately by 



AP{k) 



-8 



V93.I4 eV/ ) ^ ' 



This approximation, although not necessarily accurate on all 
scales, can be used pedagogically in the following way. This 
expression assumes fii^ ^ Qrn which is reaso nable because 
curre nt limits on neutrino energy density (e.g. lElgarav et al.l 
|2002| ) firmly exclude anything close to a Hot Dark Matter 
cosmological model. To measure the effect that neutrinos of 



a given produce we need to measure P{k) more accu- 
rately than the fractional shift given by Eqn|9] and we need 
to ensure that any systematic errors can be neglected. 

It is possible to constrain the mass of the neutrino with 
CMB experiments. However, there are strong parameter de- 
generacies that prevent this being a very precise tool. It has 
been shown, nevertheless, that even a neutrino mass as small 
as 0.05 eV could have a significant effect on the CMB tem- 
perature and po larisation power spectrum (|Hu et al.lll998l : 
lHannestadll2003l ). 

We can see in both FigO and FigO that adding two 
extra parameters {Qv, Ni,} to our fiducial cosmology yields 
different fluctuations for different values of the parameters. 
The total energy density in neutrinos £7^ is the most impor- 
tant quantity influencing P{k) (Eqn[9]l but, as explained in 
the caption to Fig. (2] the suppression of power also depends 
on A*'^ . CMB datasets are prone to serious problems with pa- 
rameter degeneracies : the bound on the neutrino mass are 
(jlchikawa et al.l |2005| ) claimed to he ^nii < 2.0 eV from 
WMAP data alone. 



2 METHODS OF MAKING FORECASTS FOR 
FUTURE COSMOLOGY EXPERIMENTS 

Cosmological experiments hold a lot of promise for helping 
to address the key unknowns in the physics of the neutrino 
sector. In this section we produce forecasts for the progress 
possible with upcoming CMB and LSS experiments. 

To make these forecasts we will use Markov Chain 
Monte Carlo (MCMC) methods to probe the n— dimensional 
parameter space of cosmological models with different neu- 
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Figure 3. The CMB temperature and polarisation power spectra for our fiducial cosmology (solid black line) compared with those with 
different values of and N^, and through Eqn[6l rrii, (coloured dashed lines). Broadly speaking, when neutrinos are added, there is 
an increase in the total mass density (at the expense of dark energy density), a change in the time of matter-radiation equality and a 
change in the amount of radiation at times before the neutrinos became cold. This produces the changes seen in this figure. There is 
a clear signal imprinted on the power spectra due to the free-streaming of neutrinos at early times (Eqn. [Sj. This is dependent on 
(Eqn. [9J and more subtly on N^, however as we will see degeneracies in the CMB analysis make it impossible to disentangle the real 
values of {A''^, Qi^} with CMB data alone. 



trino properties. Such methods (e.g. lBucher et aDl2002h pro- 
vide a powerful way of forecasting the posterior probability 
distributions for a given set of cosmological parameters to 
be probed by a future experiment. 

A MCMC is a chain of points drawn and cho- 
sen/rejected from a given distribution according to a certain 
number of rules. Points are chosen from a proposal distri- 
bution and either accepted or rejected according to the like- 
lihood attached to both the chosen point and the previous 
point only (hence a Markov chain). In the simplest algo- 
rithms, if the new point has a larger likelihood it is accepted 
and if the likelihood is lower it is accepted with a probability 
which equals the ratio of the likelihoods from the old point 
an d the new propo s al po int. The method is described fully 
by iLewis fc Bridig (|2002l ). If so, the density of points will 
describe the likelihood function in parameter space. 

Given that no real data are available, we proceed in the 
following way. We first adopt a fiducial cosmology which we 
will assume is close enough to the truth that it can be used to 



deliver simula ted data. We adopt here a theoretical frame- 
work used in iBucher et al.l (|2002h : i.e. we do not produce 
ensembles of mock datasets but assume that the simulated 
data dataaim is precisely equivalent to the prediction of the 
fiducial model with an associated, appropriately calculated, 
error bar. Any single realisation of the data would displace 
the error ellipse found but would not change its size consid- 
erably. The sum over many realisations of the data would 
produce results which are equivalent to the results obtained 
by taking the data being equal to the fiducial model, pro- 
vided, of course, the fiducial model is correct. We use this 
method for both LSS and CMB forecasts. We then attach 
theoretical errors to each fiducial data point in line with the 
errors expected for the future cosmology experiment. Then, 
we used MCMC chains to map out the n— dimensional poste- 
rior probability distribution function P{9i\datasim) , where 



P{9i\datasim) oc P{datasim\Oi) x P{9i) 



(10) 
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and P{datasim\di) is the likelihood function (given by 
exp{—x^/2) in the case of galaxy surveys where is the 
usual chi squared, i.e. the sum over all simulated data 
points of the square of the data minus the model divided 
by the error estimate at that data point) and P{9i) is the n- 
dimensional prior probability distribution function. We used 
our fiducial ACDM model with the addition of the neutrino 
parameters il^ and Ni, meaning n — 7. In the case of simu- 
lated LSS surveys more 'nuisance parameters' were needed 
to account for the way in which galaxies trace dark matter: 
two per redshift shell, making a total of 13 parameters. In 
the case of the CMB we added the optical depth of reionisa- 
tion giving a total of 8 parameters. We checked each MCMC 
chain had converged in two ways: we checked that the means 
and variances for each of the individual chains (eight MCMC 
chains were run for each case) were consistent with the other 
chains and we made a power spectrum amalysis of the chains 
l|Dunklev et al.ll2005h . Finally, we used the returned MCMC 
samples to map out the posterior probability distribution. 

Our me thods are very similar to the Fisher matrix ap- 
proach fe.g. lHu fc Haimanllioosi : ISeo fc Eisensteinll2003l ). in 
which one also assumes a fiducial cosmology and then cal- 
culates the curvature of the posterior probability, assumed 
to be an n— dimensional Gaussian, around that point. Both 
methods assume that the likelihood posterior is a slowly 
varying function of the fiducial cosmology, but our approach 
has the advantage of directly mapping out any non-Gaussian 
features of the posterior probability distribution; in the case 
of a Fisher analysis any non-Gaussianities of the posterior 
distribution function may corrupt the error analysis. 



2.1 Methods of making forecasts for future LSS 
experiments 

It has been proposed that probing the pattern of mat- 
ter fluctuations as a function of redshift yields measure- 
ments of the Hubble 'Constant' H{z) and the comov- 
ing distance D{z) via a characteristic co-moving length 
scale imprinted on the pow er spectrum P(k) via bary- 
onic oscillations (iBlake fc Glazcbrook 20031 : IHu fc HaimanI 
l2003l : ISeo fc Eisensteinll2003l ). This characteristic co-moving 
length is derivable from known physics, being essentially the 
size of the sound horizon at the time of last scattering. Given 
that matter and radiation dominate over dark energy at that 
period, its length is given by 



Jo 



1 



da 



(3(l + i?))i/2 aE{a) 



(11) 



where R — 3pi,/4p^ introduces a weak dependence of the 
sound speed of the plasma on Q,b and Srec is the scale factor 
at the redshift of recombination. We can interpret this as a 
sound speed that is equal to c/\/3 in the case of radiation 
only and which decreases as a larger fraction of baryons is 
present. 

Therefore we know the precise length of this standard 
rod for a given cosmology. Furthermore, this distance is de- 
pendent on h, Qrn and weakly on Qb but it is very insensitive 
to w given that this acoustic scale was set up in the early 
universe when it is thought that dark energy plays a small 
role. 



A related geomet ric probe of the evolution of the Uni- 
verse was proposed bv lAlcock fc Paczvnskil (|l979l ) by assum- 
ing that any characteristic length present in the Universe, 
and which can be seen in both radial and angular direc- 
tions in the sky, produces a measurement of the product 
H{z)Da{z) [where Da{z) = D[z)/{1 + z) is the angular di- 
ameter distance] simply because, by isotropy, the length has 
to be the same in both directions. 

This proposed wigg l es test (e.g 'Blake fc Glazebrookl 
l2003l : lHu fc Haimanl200a : ISeo fc Eisenstcin 2003i ') can there- 
fore simply be considered as a more powerful version of 
the Alcock-Paczynski (AP) test and should be more use- 
ful in constraining dark energy as a geometrical test which 
is weakly dependent on w. 

More precisely, the comoving sizes of any object or any 
feature has a transverse {r±) and parallel ry projection on 
the sky which can be used for this kind of test. These comov- 
ing features relate to the sizes of the observed angular and 
redshift distances A6 and Az respectively via the following 
relations 



cAz 

- {i + z)DAiz)Ae. 



(12) 



Hence, when the true scales are known the measure- 
ment of the pattern in z or the line of sight and 9 or in 
the plane perpendicular to the line of sight yield a measure- 
ment of H{z) and Da{z), whereas if the precise value of r 
is not known, it is only the product H{z)Da{z) that can be 
measured. 

To be able to estimate the errors on cosmological pa- 
rameters we have to know how well we will be able to mea- 
sure the power spectrum. Under Gaussian conditions, the 
fractional error on the power spectrum will depend on the 
total volume of the survey as this will determine the num- 
ber of Fourier modes that will be sampled by the survey; the 
associated error is usually called cosmic variance. Further- 
more there will be an uncertainty given by the finite number 
of galaxies. This is commonly termed shot noise. Hence the 
fractional error on t he power spectrum can be written as 
(|Feldman et al.ll 19941 ) 



47r 



\ P J Vk^AkAfj. 



l + nP 
nP 



(13) 



where P is the power spectrum measured, n is the comoving 
number density of galaxies probed and fi is the cosine of the 
angle defined by our line of sight and the line joining pairs 
of galaxies in 3D space. 

If we make a measurement of P{k) several effects can 
change, not only the shape, but also the height of the power 
spectrum, as a function of redshift. If a wrong cosmology is 
assumed our estimates of distances and volumes are incor- 
rect. This produces a distortion that is visible in the power 
spectrum as rings in the (perpendicular, parallel) k plane 
(|Hu fc Haimanll2003h . Furthermore, at a higher redshift, the 
growth factor changes the height of the power spectrum. 
However if a wrong cosmology is assumed, then we have 
wrong measurements of the cosmological volume in which 
we made our survey and therefore the height of the power 
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spectrum is measured incorrectly. These distortion effe cts 
have been thoroughly described in iBaUinger et~aLl (|f996l l. 

If the galaxies are measured in redshift and we do not 
have any information about their peculiar velocities, we then 
retrieve a redshift space power spectrum. This will differ 
from the real space power spectrum in two ways, firstly 
we obtain more correlations at low k because large-scale 
bulk flows point directly towards matter overdensities due to 
gravitational pull. The refore there is an enhancement factor 
(1 + fkfi'^f derived in lKaiseij (|l987l ). where ^ is the cosine 
of the angle that the mode which is probed makes with the 
line of sight, and fk is the derivative of the natural log of 
the over-density with respect to the natural log of the scale 
factor /fc = dln5/dlna. At large scales, structures appear to 
be closer to each other in redshift space. On the other hand 
at small scales, given the circular velocities of satellite and 
relaxed structures, redshift-space structures appear to be 
elongated along the line-of-sight. This c reated the so called 
finger-of-god effect jPeacock et al.|[200ll ). Here we are inter- 
ested in the large-scale fluctuations, and we simply ignore 
the small-scale power , i.e. we include the first effect but not 
the finger-of-god effect. Note that a complete analysis with 
real data will need to take both effects into account. 

When a galaxy survey is performed it is vital to distin- 
guish between the correlations that we measure in the galaxy 
distribution and how that distribution relates to the distri- 
bution of dark matter. We assume throughout that the non- 
linear galaxy power spectrum is related to the linear matter 
power spectrum via Pgik) = h^Pjk) + Pshot- This is moti- 
vated by the halo model (e.g. ISeliakllioOOl : iPeacock fc SmithI 
I2OOOI ). a galaxy population will most likely not trace matter 
in an unbiased way and given their discreteness and number 
density there will be an extra shot noise power attached to 
the power spectrum. 

Therefore in order to estimate cosmological parameters 
from a galaxy catalogue we first choose a fiducial cosmol- 
ogy which will allow us to change the coordinates from 
(RA,DEC,z) to spatial coordinates in redshift space. If, pre- 
sumably by luck, we choose the real cosmology that governs 
our Universe, then the power spectrum that we will measure 
as a result of the galaxy autocorrelation function is 



Pobs{k, fi) 



1 + 



P{k) + Pshot 



(14) 



However if we chose an incorrect reference cosmology, 
then we will observe the signal produced from the power 
spectrum derived from the real cosmological parameters dis- 
torted by the incorrect assumption we have made. This re- 
sults in the following expression, for what measurements of 
P{k) we would retrieve 



order to do this, with real data, we simply need to repeat 
this process for may different cosmological parameters and 
find which cosmological parameters do not produce any in- 
consistency by comparing predicted models and measured 
models assuming the cosmology we are testing. 

In both equations [T?] and [T5] a relationship between 
(fc,/i) and (fcrcf , Mrcf ) depends on cosmology in the follow- 
ing way. 



1.2 _ i.2 I 7.2 J.2 _ 7,2 I 1.2 

f- — fi-x ^ '''II 'S'ref " '^±,ref ""ILrof 

k\\ — Ilk fc||_ref = Mreffcrof 



(16) 



which produce the following relations that we use to convert 
from a given angle and k value in the sky to another dis- 
torted angle and k value if the incorrect cosmology is chosen 
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2 

Mrcf 



+ (1 

Hiz) 



/ircf 



PA.rof(^)gA.rof(-) 



(17) 



This is a thin shell approximation and strictly speaking 
should not be used in this analysis. However, we argue that 
the effect on the power spectrum of using a thin shell ap- 
proximation, if the shells used are thick, is equivalent to 
considering a power spectrum convolved with a narrow k 
space window function. We therefore consider, for forecast- 
in g purposes, this is a good en ough approximation as shown 
bv lciazebrook fc Blakel (|2005l ) for similar survey geometries. 

Here we will produce forecasts using Markov Chain 
Monte Carlo (MCMC) methods to predict the posterior 
probability distribution for a given set of cosmological pa- 
rameters. In this case, given that no real data are available, 
we proceed in the following way. We chose a fiducial cos- 
mology which we will assume to be close to what we think 
reproduces the real distribution of galaxies. In order to pro- 
duce a MCMC that will be able to provide us with the er- 
ror forecasts we need to have a prescription to compare the 
log(likelihood) difference between two potential models. We 
compare a model with another model which has been dis- 
torted by the incorrect choice of cosmology. In other words 
in order to produce a log(likelihood) difference between two 
models we compare the results from Eqn[T4] with results 
from Eqn[l5] having assumed that the errors in the power 
spectrum will be given by Eqnll3l We calculate the for 
simulated galaxy surveys by using a method involving bins 
in fi, redshift z and k, and including only k modes up to a 
conservative value of kmax = 0.2 h Mpc~^, beyond which 
the power spectrum becomes non-linear. 



Pobs {k^ci , /^rcf ) 



D\{z),,iH{z) 2 

~ '-'true 



D\{z)H{z) 



1 + 



rue 

{z) Pfcncik) + -Pshot- 



(15) 



Hence if an incorrect cosmology is assumed there will be 
an inconsistency when it comes to compare the power spec- 
trum expected given the assumptions we made about the 
cosmological parameters. With the power spectrum derived 
from the data given assumed distances, we will be able to 
use this inconsistency to probe cosmological parameters. In 



2.2 Methods for CMB (Planck) datasets 

We will briefly describe how we use MCMC meth- 
ods to produce forecasts for future CMB experiments 
and specifically for the angular power spectrum mea- 
surement to be achieved by the CMB satellite Planck 
( http: / / www.rssd.esa.int/Planck ) . Planck, like WMAP be- 
fore it ijSpergel et al.|[200a . I2OO6I ') . will produce a CMB map 
by using several radio bands to remove foregrounds, espe- 
cially dust and synchrotron emission arising from the Milky 



© 0000 RAS, MNRAS 000, 000-000 



Determining neutrino properties using future galaxy redshift surveys 7 



Way and point sources. The nearly constant background cor- 
responding to a 2.7 K black body will be removed, as will 
the pure dipole component due to motion with respect to 
the CMB frame. This leaves a map of CMB anisotropies 
Q{6,(j>) which is decomposed into spherical harmonics 



(l.m) 



(18) 



The statistical properties of the coefficients aim are 
translated into constraints on the underlying cosmological 
parameters which produce the fluctuations and determine 
their evolution. Theories for cosmic inflation (e.g. iGuthI 
Il98ll ) predict that the primordial density fields are Gaus- 
sian in nature, in which case the average (flim) = and it is 
the variance 



{|ciimCl*/„j/ I) = SlmSl'm'Cl 



(19) 



that contains the crucial cosmological information. 

The values of Ci can be calculated fro m a 'Boltz- 
mann code' such as CAMB (|Lewis et al.ll2000l ') or CMBfast 
l|Seliak fc ZaldarriagalligQel ) ; in this work we use CAMB. 

To produce forecasts, we need to have a good idea of the 
errors. Even in the absence of any experimental errors, there 
is an intrinsic 'cosmic variance' because there is a limited 
number of modes (i.e. a finite number of aim for a given /) 
that can be measured on the sky. The total number of aim 
is simply related to the total number of independent pixels 
available in the sky anisotropy map. The cosmic variance 
error can be written as 



AC; 

Ci 



(20) 



where the factor 21 + 1 is the number of modes aim measur- 
able from the data which is proportional to the fractional 
area of sky /sky studied and the factor 2 reflects that fact 
that the direction of the modes on the sky is unimportant. 

We must also account for experimental sources of error 
on the aim measurements. If we assume that the sources of 
variance are independent, which is obviously fine since one is 
cosmological and the other experimental, then the variances 
should be added in quadrature to get the total error estimate 



AC, 



fsky{2l + 1 



■{Ci + Nf) 



(21) 



where Ni will be defined as a function of the experiment. 

The experimental error will be determined by the sen- 
sitivity CTc and the angular resolution (beam size) Of,. There- 
fore, if we have a certain number of frequency channels ob- 
serving the CMB for a given experiment, the noise will be 
given by the following expression (,Eiscnstcin ct al.nlQQQl ') 



N, 



channels 



i{i + i)et, 

81n2 



(22) 



assuming that the beam smears the measurement on a given 
/ scale. 



The sensitivity of the experiment will depend on the 
experimental setup through the equations 



NET 



sky 



(23) 



where NET is the noise effective temperature in one second 
for each of Udet detectors, t is the integration time and 6sky 
is the solid angle of the survey . 

The probability of a value of aim given a Ci value can 
now be written as 



p{aim\Ci) 



V2^i 



:exp 



2 

2Ci 



(24) 



so the probability that a given realisation of the sky is pro- 
duced, given a set of cosmological parameters, is simply the 
product of all the aim probabilities which, in the case where 
the values of aim are not correlated, gives 



pisky\C'i) = TT -— ^ exp 
,™ V2tvCi 



2 

2Ci 



(25) 



We can then use Bayes theorem to find that 
jBucher et al.ll2002l ') 



l°g(^) -log 



p{aim\A) 
p{aim\B) 



log 



Ci,A + N^ 
C,,B + Nf 

Cl.A + 



+ 



(26) 



Ci,B+Nf 

However, in the case of future CMB observations like 
those from Planck, we will have access to data which 
gives information on both temperature and polarisation 
anisotropies. It is common to separate the polarisation 
modes in the CMB as E and B modes which are related 
to 'curl free' and 'curl' vector fields. This separation is done 
because most primordial effects generate only E modes, and 
B modes are only generated if there is a gravitational wave 
background or by gravitational lensing of E modes. Much 
additional cosmological information can be extracted from 
the E-mode power spectrum as well as the cross correlation 
between the E and T (total power) spectra. 

Mathematically, Eqn[26]is modified in the following way 
in the case where the model describes all the power spectra 
produced from rand om Gaussian fields w ith a certain degree 
of cross correlation (|Bucher et al.|[20o3 ). 



log 



Pa\ 
Pb) 



log 



p{aim\A) 
p{aim\B) 



) + 



= ^ I](2« + 1) {Tr{I - MaMb' 

log [DetiMAMg^))) 



(27) 



where / is the identity matrix and the matrices AIa and A/ b 
are given by the values of the individual power spectra and 
their cross correlations at a given mode / following 



M 



Cr + Ni 



rprp2 



(jTE 

cr + Ni' 



EE^ 



(28) 
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Figure 4. Posterior probability distribution from an MCMC run determining the quality of measurements that the Planck CMB mission 
will have in measuring the absolute neutrino mass scale. Contours are one sigma and two sigma levels obtained by taking the 68% and 
95% of the samples which best fit the fiducial model. It is clear that degeneracies are unlikely to allow us to get a clear signal arising 
from a massive neutrino from CMB data alone. Data from Planck will improve current constraints on the neutrino mass even though a 
direct detection seems unlikely as many models with = are not rejected at the one sigma level. 



We can therefore produce a forecast for any CMB ex- 
periment provided we have access to rehable predictions of 
the noise level of the experiment for both temperature and 
polarisation. Here, we produce forecasts involving the entire 
likelihood for a Planck-like experiment. The Planck specifi- 
cations used here are given in Table[l] 

We plot in Fig|4]the results we find for an MCMC run 
assuming Planck data alone. As we can see the main degen- 
eracies will remain in Planck data, noting the strong nega- 
tive correlation of Q^, and as which can be explained by less 
growth with a cosmology with neutrinos. Specially models in 
which f2y = 0.0 cannot be ruled out if values of m„ turn out 
to be near the bottom range allowed by particle physics. The 
upper limits derived from Planck data alone are better than 



Frequency 


70 


100 


143 


217 


Beam width 0f, / arcsec 


14.0 


9.5 


7.1 


5.0 


NET /(fiKy^) 


212(300) 


56(80) 


56(80) 


84(120) 


Detector number n^jej 


12 


8 


12(8) 


12(8) 



Table 1. Experimental limits used in the forecasts of future CMB 
(Planck) data. Values in parenthesis are for polarisation exper- 
iments. We have assumed one year integration time and a sky 
coverage of f^ky = 0.65. These data are taken from the Planck 
blue book ( http: / /www.rssd.esa.int/Planckj l 

current estimates for the neutrino mass: marginalising over 
other parameters they lie around Am^ ~ 0.65 eV (1 a). 
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This result is roughl y in accordance with other forecasts 
(|Hu fc Tegmarkll 19991 ). where the authors find a lower value 
of Am^ ~ 0.25 eV (1 a) using a Fisher matrix a n alysis and 
broadly similar assumptions. In IHu fc TegmarkI (|l999t ) the 
assumptions of noise of the Planck experiment were slightly 
different than the ones assumed here, and their Fisher anal- 
ysis does not a ccount of any non-Ga ussian distribution of 
the parameters. Ilchikawa et al.l (|2005h claim that if neutri- 
nos have a mass scale smaller than ~ 1 eV then there is very 
little useful constraining information in CMB datasets. 

We point out here that we have made an important 
approximation in our analysis. We have assumed that we 
are in the limit where we have such a small mass for the 
lowest mass eigenstate, that we can set the masses of the 
lowest mass eigenstates to zero (the two lowest for a normal 
hierarchy, and the lowest only for an inverted hierarchy; see 
Fig. [1} . This approximation is valid if the masses are small 
enough that we are not in the quasi-degenerate scenario. 
Given that all the scenarios we consider here assume a mass 
below 0.25 eV, this is a reasonable assumption. We have not 
considered here the possibility that different neutrino masses 
may be measurable by cosmological data as this is beyond 
the capabilities of the surveys we consider. 



2.3 Results from LSS and effects of priors. 

It was pointed out in lElgarOv fc Lahavl (|2003l ) that the role 
of priors is vital in retrieving information about the neutrino 
mass from cosmological surveys. This is because of the large 
set of parameters that determines the behaviour of cosmo- 
logical perturbations and the complicated relation between 
them. There is for instance a tight correlation between the 
matter density and the neutrino parameter fi^. This occurs 
because over most of the evolution of the Universe, neutri- 
nos behave just as cold dark matter, hence it is necessary 
to know the cold dark matter density well to determine ef- 
fectively by how much it effectively changes as neutrinos 
become non-relativistic. 

Given that one of the major effects of a high Qi, is a 
strong damping of the power at small scales due to free 
streaming, it is strongly correlated with the value of n^. It is 
very important to measure Us in order to deduce a neutrino 
mass. It is possible that the discrepancies between current 
results that do agree broadly but not in the detail, comes 
down to the priors used, for instance (|Elgar0v et al.ll2002l ) 
did not assume a varying scalar index in their analysis. They 
chose to give results for different choices of = 0.9, 1.0, 1.1. 
One may use theoretical priors to consider only values close 
to ns = 1 but in doing so one must be wary of the large 
degeneracy between Us and ^Iv 

Most current results rely on priors in one way or an- 
other. We argue that in order to obtain a clean measurement 
of the neutrino mass it will be important to use cosmologi- 
cal data in a way that we are sure that our result in not too 
prior dependent. We will argue in this paper that a future 
CMB data set, such as one from Planck, combined with a 
galaxy survey, such as one from the SKA, can provide the 
data necessary for a robust analysis of neutrinos via cosmol- 
ogy. 




k / h Mpc"^ 



Figure 5. Illustration of the sensitivity a future LSS (SKA) sur- 
vey will need to measure and constrain the absolute neutrino 
mass scale m,,. The error bars illustrate the accuracy of the LSS 
measurement assuming the validity of Ean ll3l with P{k) taken 
from a fiducial cosmological model which neglects massive neu- 
trinos {Ut, flc, w, h, Us, (T8}= {0.04, 0.26, -1, 0.72, 1.0, 0.9}. 
Note that the plotted error bars should be totally uncorrclated 
because of the large area coverage and redshift depth, and hence 
broad real-space window function, of the SKA survey. The solid, 
dashed and dotted lines correspond to the addition to the fiducial 
model of neutrinos of mass 0.05 eV, 0.25 eV and 0.5 eV (taking 
Ny = 3), all of which are still possible values given the constraints 
from current data sets. The diff'erent colours represent indepen- 
dent SKA measurements in three redshift slices: 0.5 2 < 1.0 
(blue), 1.0 ^ 2 < 1.5 (green) and 1.5 ^ 2 < 2 (red). 



3 MEASURING NEUTRINO MASSES WITH 
THE SKA 

3.1 A back-of-the-envelope calculation 

As we have outlined, a measurement of the shape of the 
power spectrum at large scales as given by Eqn|9] would be 
able to constrain for the neutrino sector. Here we at- 
tempt to estimate the extent to which we would be able 
to measure this effect with a future large-scale-structure 
(LSS) survey. In this section we outline a simple back-of-the- 
envelope calculation to estimate the cosmic volume required, 
and then in Sec l3.3l we undertake a more detailed calcula- 
tion that takes into account the full effect that neutrinos 
have both on LSS and CMB data. 

However, data produced by redshift surveys on large 
scales is often cosmic- variance limited. In this case, the frac- 
tional error with which we can measure P{k) is proportional 
to the number of k modes present in the survey volume, and 
hence the accuracy is inversely proportional to the square 
root of the cosmic volume fEqn ll3|l . 

If we ignore for now the role of priors that, in real- 
ity, complicate the analysis (Sec. 13.3) . the upper limits from 
surveys such as the 2dF Galaxy Redshift Survey (2dFGRS) 
and the SL OAN Digital Sky Sur vey (SDSS) are around 
nil, ~ leV (jTegmark et al.l |2006| ) . Hence, a survey with 
~ 400 times the cosmic volume would have errors on large 
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Figure 6. Analysis of tlic accuracy that a future SKA survey will have in constraining a neutrino mass scale. Contours arc one sigma 
and two sigma levels obtained by taking the 68% and 95% of the samples which best fit the fiducial model. The SKA alone will be able 
to, by accurately measuring the shape of the power spectrum in several different redshift bins, measure a signal from the neutrino energy 
density at more than the three sigma level by marginalising over other parameters and without any further priors needed by current 
data sets. 



scales a factor of ~ 20-times smaller and therefore would 
be able to probe ~ 0.05 eV at the bottom end of the 
range allowed by particle physics, provided the relation given 
in Eqn[5]holds. Such an experiment would, by detecting the 
signal from neutrinos imprinted on the large-scale structure, 
be able to probe the entire neutrino sector allowed by cur- 
rent particle physics experiments 

Furthermore, if we do have such a large volume to search 
for features imprinted on the galaxy power spectrum, other 
cosmological parameters will also be accurately determined. 
The problems with priors and parameter degeneracies that 
plague current analyses should be much reduced. Current 
redshift surveys also suffer from the problem of correlated 
errors on large scales because window functions of these sur- 



veys are necessarily narrow in real space, and hence broad in 
k space. An 'all sky' survey reaching to very high redshifts 
would be the optimal way to probe such signals, but we will 
consider here an, as yet hypothetical, 20, 000 deg^ LSS sur- 
vey reaching to redshift z ~ 2. Such a survey would provide 
the huge increase in cosmic volume needed to comprehen- 
sively probe neutrino properties. For specific calculations we 
will use simulated LSS datasets from surveys with the SKA, 
although clearly our conclusions will be valid for any LSS 
survey with similar reach in sky area, redshift depth and 
galaxy number density. Further details concerning future 
redshift surveys with the SKA can be found in the foUow- 
in g references: details concer ning galaxy number densities 
in lAbdalla fc Rawlingsl (|2005h : details concerning P{k) mea- 
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Figure 7. Analysis of the accuracy that a future SKA survey will 
have in constraining an absolute neutrino mass ~ 0.25 eV. 
The plot shows the forecast of the unnormalised probability one 
expects to have for the sum of the neutrino masses. Here the fidu- 
cial value has been chosen to be Em^ = 0.25 eV. The vertical red 
dashed lines correspond to one, two and three sigma confidence 
levels. By accurately measuring the shape of the power spectrum 
in three different redshift bins, the SKA alone will be able to, mea- 
sure a signal from the neutrino mass at more than three sigma 
level by marginalising over other parameters and without any 
further priors as are needed for current data sets. 



surement in lAbdalla et al] (|2006|): a comparison w ith other 
future LSS datasets in lRawlings fc Abdallal l|2006h . 



3.2 



The cosmological imprint of neutrinos in a 
future LSS data set 



We illustrate in Fig[S] the promise of a future LSS (SKA) 
survey for detecting the imprint of neutrinos on the power 
spectrum P{k). We note that we have chosen bins in red- 
shift space which correspond to intervals in comoving space 
larger than the A corresponding to fc = O.Olft^^Mpc, hence 
the theoretical error bars plotted are uncorrelated. Fig[5]is 
the result of a simple analysis in which all other cosmo- 
logical parameters are fixed at fiducial values and is 
varied through the neutrino density parameter Sl^ (using 
Eqn|6] in a quasi-degenerate scenario with = 3). We 
can see clearly that this agrees well with the back-of-the- 
envelope calculation presented in Sec l3.1l with even models 
with ~ 0.05 eV starting to show significant deviations 
from the fiducial model which ignores neutrinos. 

Another key point illustrated by Fig[5]is the importance 
of undertaking a LSS survey in independent redshift shells. 
A common criticism of attempts to constrain rrii, from LSS 
measurements is that the feature being looked for maybe 
being masked, or even mimicked, by features introduced by 
the use of galaxies as an indirect probe of the underlying 
fiuctuations. A LSS survey which can probe several redshift 
shells can give a very good handle on any such systematic 
errors present in a galaxy redshift survey. One specific worry 
that has to be addressed is whether any biasing effect from 



the galaxy power spectrum would be able to produce any 
scale-dependent mo dulation mim icking the damping tail of 
massive neutrinos (|Seliak|[2000l ). First, any worry on this 
specific issue is attenuate d from the t h eoretical expecta- 
tion from the halo model dSeliakl (I2OO0I ): IPeacock fc SmithI 
(|2000l ): and see lAbdalla et al.l (|2006l ))~that the galaxy bias 
is, at least on large scales, reliably modelled as a constant 
multiplicative factor plus a constant additive shot noise. It is 
only at small scales that this assumption should break down. 
Hence, the neutrino mass estimate should be robust against 
this kind of systematic. Second, as illustrated in Fig[5] the 
huge number of galaxies of an SKA-like LSS survey would 
allow us to probe the galaxy power spectrum with several 
galaxy types in each redshift shell, and therefore produce 
an independent power spectrum for each galaxy type. If we 
are concerned about any systematic problem due to galaxy 
bias, we should be able to compare the power spectrum for 
different galaxy types in each shell. All shells should have 
the same signal arising from neutrino physics, whereas each 
power spectrum would be different if there is an effect mim- 
icking the effects of neutrinos because each is made up of 
samples of galaxies wh i ch ha ve different bias and clustering 
properties. Cole et al] (|2006l ) have shown that it is impor- 
tant to correctly model any scale dependent bias which is 
luminosity dependent. With a future LSS survey it will be 
important to have enough galaxies to correctly model this 
effect. With an SKA survey the number of galaxies will be 
large enough so that we are able to model the scale depen- 
dence of the bias as a function of galaxy properties which 
we can retrieve from the data, i.e. as a function of hydro- 
gen mass and circular velocity which should correlate to the 
dark matter mass of the galaxy. There is therefore a bet- 
ter prospect to model a scale dependent bias than may be 
possible with optical surveys. 

However, there is another potential problem that Fig[5] 
fails to address. We still need to test whether constraints 
on neutrino parameters are robust to problems induced by 
degeneracies between cosmological and neutrino parame- 
ters. Such degeneracies could make an LSS experiment more 
heavily reliant on priors than on the data itself, and hence 
yield misleading forecasts. 



3.3 MCMC methods probing a quasi-degenerate 
scenario with the SKA. 

In order to probe the cosmological parameter space for these 
models we use a Markov Chain Monte Carlo (MCMC) meth- 
ods in order to obtain a prediction of the posterior proba- 
bility for a given model and a given future survey. We use 
a AC DM model with neutrinos as a model cosmology and 
use the following 13 parameters in our MCMC chains: {Q,b, 
Q.C, fii/, h, Us, erg, A^^, 61, pi, 62, P2, &3,P3}. The first seven 
parameters are the usual cosmic parameters and the last six 
parameters are the multiplicative bias and an additive shot 
noise power Pshot arising from the halo model, in each one 
of the three redshift bins, that we marginalise over. 

We will focus on two MCMC-based studies: one which 
assumes the absolute mass scale for the neutrinos ~ 
0.25 eV chosen to be the lowest mass which an SKA survey 
can, on its own, measure sufficiently to discriminate between 
a normal and inverted hierarchy; another which combines 
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Figure 8. We illustrate the main degeneracies present in LSS 
(SKA) and CMB (Planck) data and how the combination of both 
approaches can improve the errors on Q^. We plot MCMC sam- 
ples from both experiments, separately simulated, with different 
parameters along the x and y axis. The top panels illustrate that 
the CMB will measure Us accurately whereas there is a degen- 
eracy between and ris in LSS data. The bottom panel shows 
that there is an anti-correlation between the parameters erg and 
Qv for CMB data whereas there is a positive correlation in the 
case of LSS data; this is due to different effects: for the CMB, the 
Qjy - (Tg degeneracy is probing the amount of growth from recom- 
bination to the present day; for the LSS, the Qi, - eg degeneracy 
is due to the shape of the power spectrum which is determined 
by Oi/ and which is central to the determination of erg. 



Planck and SKA data on the assumption that the mass scale 
is nil, = 0.05 eV, the lowest allowed by particle physics. 

3.3.1 SKA-only study at Enn ~ 0.25 eV. 

We first choose to examine whether a survey with the SKA 
will be able to probe as low as ~ 0.25 eV and be able to 



completely rule out models which have a quasi-degenerate 
neutrino mass spectrum. We therefore choose a fiducial 
model which has only two massive neutrinos with the fol- 
lowing cosmological parameters {fib, He, ^i^, h, Us, as, N,^} 
= {0.04,0.255,0.005,0.72,1.0,0.9,2.0}. We impose an uniform 
top-hat prior on < N^, < 3 assuming that the standard 
model for neutrinos imposes three neutrino eigenstates. The 
fiducial bias for all the redshift bins has been conservatively 
set to one but we note that if a larger bias is found, which 
is likely to be the case, especially in higher redshift bins the 
significance of the measurements of neutrinos will improve. 
We assume a fiducial shot noise value equal to one over the 
number density of galaxies. 

We plot in FiglH] and FigE] the MCMC samples for an 
analysis following the methods outlined. We note that with 
this future SKA survey alone, there should be a ^ So" detec- 
tion of the absolute mass scale of the neutrino. The measure- 
ment of the power spectrum is accurate enough that other 
parameters are well constrained so stringent priors are there- 
fore unnecessary for a detection of the neutrino signature to 
be made. We argue that this is of extreme importance for a 
result to be considered robust not only by the astrophysical 
community but also by the particle physics community. For 
such a sample we can estimate the direct forecasted poste- 
rior on the sum of the neutrino masses. We plot this in FiglT] 
from which the neutrino absolute mass scale can be detected 
at more than the three sigma level. 

3.3.2 SKA plus CMB study atY^rm^ 0.25 eV. 

It is important to understand and predict how other data 
will affect the predictions of Sec l3.3.l1 Even though we have 
argued that restrictive priors may be unhelpful when it 
comes to interpreting how much information is encoded in 
a given data set, it is certain that additional data sets can 
provide a huge gain in sensitivity for a given cosmological 
experiment as has been pointed out by many authors (e.g. 
IHu et al.|[l993 ). 

If we ignore the problems of parameter degeneracy 
a LSS survey which can get redshifts for sources over 
20, 000 deg"^ out to a redshift of two would be able to reach 
rriv ~ 0.05 eV (Fig[2J. However the problem of parameter 
degeneracy plagues any interpretation of these data. For in- 
stance if we are assessing models which have small values of 
fi^ then the suppression of the power at small scales is very 
small and can be mimicked by a change in Us in the primor- 
dial power spectrum. Hence a LSS survey would be unable 
to distinguish between a signal left over from an infiationary 
phase of the Universe and a signal due to neutrinos. 

In Fig|9] we see the full posterior probability distribu- 
tions for the combination of LSS (SKA) and CMB (Planck) 
experiments. We illustrate in Fig|8]the Q.^ - Us degeneracy 
present in a future SKA survey on its own. This degeneracy 
is made worst because we have included k modes up to a 
conservative value of kmax of 0.2 h Mpc~^. The power spec- 
trum becomes non-linear at larger scales (smaller k) than 
kmax and if we could find a satisfactory way of modelling 
these non-linearities we would be able to use these data to 
help break this degeneracy. However Planck will be able to 
probe high k values by measuring high I values in the CMB 
power spectrum. As gravitational growth is much less ad- 
vanced ai z = 1000, non-linearities are not an issue and 
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Figure 9. Similar figure to Fig|6]for a LSS (SKA) plus CMB (Planck) experiment. Contours are one sigma and two sigma levels obtained 
by taking the 68% and 95% of the samples which best fit the fiducial model. The combination of these two experiments will be able to 
constrain several cosmological parameters to the 1 per cent level. Most notably, the energy density in neutrinos will be constrained to 
~ 0.015 eV and the number of massive neutrino species to ~ 0.5. 



CMB data are therefore also ideal to constrain Us- As we 
can see from Fig|8] Planck will measure very well inde- 
pendently of most other parameters. Hence a combination 
of these two data sets will improve the error bars on the neu- 
trino mass by a factor of five assuming a fiducial model with 
= 0.25 eV. We plot the improvement on the estimates 
in FigfUJI 

Another degeneracy that plagues estimates of the neu- 
trino mass via cosmological methods is the value of the ini- 
tial fluctuations which can be parameterised as the value 
it has at the CMB via As or the current strength of the 
anisotropies parametrised by erg. Given that one of the ef- 
fects of a massive neutrino in cosmology is the change in the 
growth factor for a lengthy period of the history of the Uni- 
verse, a direct comparison between the scale of anisotropies 



on the CMB and more locally with redshift surveys can also 
improve greatly the estimates of the neutrino mass. 

We plot in Fig[S]the as - ^1^, degeneracy for CMB and 
LSS surveys. Whereas for the CMB, growth is responsible for 
the large degeneracy between Q^, and as, for the a LSS sur- 
vey is is the shape of the power spectrum which is mainly re- 
sponsible for the observed degeneracy, i.e. erg fixes a weighted 
integral under P{k), whereas Q^, changes the overall shape of 
P(k). By combining the two estimates a huge improvement 
can be achieved. 

However a poor knowledge of the bias parameter would 
not allow us to combine these data sets in a consistent way 
and the improvement would be much smaller. We have as- 
sumed that the bias of high redshift galaxies can be mea- 
sured and marginalised via measurements of modes with 
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Figure 10. Confidence levels obtained for the neutrino mass as- 
suming mi, = 0.25 eV for LSS (SKA) and CMB (Planck) data. A 
comparison between this figure and Fig[7]shows how the accuracy 
of the measurement is improved by a factor of five when both ex- 
periments are considered. Vertical red dashed lines correspond to 
one, two and three sigma confidence levels. 



different sky orientation, through the influence of redshift 
space distortions (see Ean ll4|l . In our MCMC chains the 
bias is constrained to a few per cent level and marginalised 
over. For a more realistic simulation we could include other 
effects such as small scale redshift space distortions. Current 
spectroscopic surveys with much smaller samples are able to 
measure the bias to about the several per cent level by mea- 
surin g bias using higher order statistics such as the bispec- 
trum lVerde et all (|2002D . so we argue that the assumptions 
we have made here are not over-optimistic. 

However we point out that it is possible that the data 
at high redshift will allow us to probe linear scales up to a 
much larger k values in which case the constraints coming 
from LSS alone would be greatly improved. Also a better 
understanding of non-linearitie s coming from the ha lo model 
as well as N-body simulations (|Springel et al.|[2005h are very 
promising, and there are realistic hopes that we could use 
data from the mildly non-linear regime without introducing 
large systematic errors and hence improve the constraints 
shown here. 

3.3.3 Directly probing neutrino hierarchies with 
cosmological data with "^nii ^ 0.25 eV. 

In our analysis of mock future data we have assumed that 
the number of massive neutrinos Ni, is constant. If we as- 
sume the standard scenario for neutrinos this can arise in 
a hierarchical case when the hierarchy is either an inverse 
hierarchy or a normal hierarchy (FigdJ. If the hierarchy is 
quasi degenerate, then the number of massive neutrinos will 
be characterised by A'^^ = 3 and this represent a density of 
neutrinos equal to 422 million particles per m^. 

As we stated in the Sec l3.3.2l we have assumed in our 
fiducial model a number of neutrino A^^ = 2 which corre- 
sponds to an inverted hierarchy where the neutrinos are not 
all degenerate in mass (Fig[T]). By assuming that is a 




Figure 11. Confidence limits on neutrino hierarchies with LSS 
(SKA) and CMB (Planck) data. Wc assess whether a future sur- 
vey can measure the number density of massive neutrinos and 
hence have a direct measurement of the hierarchy of the neu- 
trino mass scales assuming a total energy density in relativistic 
particles to be = 0.005 which corresponds roughly to a neu- 
trino mass of 0.25 eV. The vertical red dashed lines correspond 
to one, two and three sigma confidence levels. We can see that 
the parameter A^ is a difficult parameter to measure and oven 
with exquisite data coming from future galaxy surveys and future 
CMB experiments it will be hard to constrain this parameter to 
better than a &Ny > 1 . For this specific scenario the value Ni, = 1 
can be rejected at 2.7 a. 



continuous variable we can therefore compare and contrast 
models which assume a normal hierarchy and models that 
have an inverted hierarchy. A model with — 1 would 
correspond to an normal hierarchy. 

As we ran MCMC chains for LSS (SKA) data alone and 
CMB (Planck) data alone we noted that for a top-hat prior 
< A^i^ < 3, which is what one would expect in a standard 
neutrino scenario without sterile neutrinos or any other ex- 
otic relativistic particle that could contribute significantly 
to the energy density of the Universe, both of these exper- 
iments do not constrain A''^ significantly. In fact as we can 
see from the relevant panels of Fig[6] and Fig|4j the prior 
plays a bigger role in the posterior than the experimental 
information itself if we consider each experiment indepen- 
dently. It would be more instructive if this were real data to 
widen the prior allowing more exotic scenarios to be allowed 
by the model. However, when we combine both experiments 
we obtain a constraint which is slightly prior dependent but 
which is now mostly determined by data and this combina- 
tion of data would yield a measurement of A^^ — 2±0.5 (la) 
and would reject the value Ni, = 1 at just under 3a. We plot 
in Fig llll the expected posterior probability for the value of 
Ni, in this case. 

This combination of and Ni, would yield invaluable 
information in order for us to be able to disentangle the 
origin of the neutrino mass and have a convincing theory 
which would explain the huge difference between the masses 
of other particles and the neutrino mass. We plot in Fig |12l 
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Figure 12. Confidence limits on neutrino hierarchies with LSS 
(SKA) and CMB (Planck) data. Illustration of the distribution of 
MCMC samples in the {Af^, f2^} plane. The covariance around 
the fiducial {2, 0.005} point is such that it is possible for us to 
discriminate between a normal and inverted hierarchy. For higher 
values of a fiducial Qi, the detection is clearer, and as we can see 
from Fie llSI for lower values of a fiducial Qi, the error becomes 
larger. 




Figure 15. Wc show here MCMC samples in the 2D plane {N,j, 
fit/}: the fiducial model chosen has a neutrino mass of 0.05 eV and 
2 massive neutrino cigenstates, and both CMB and LSS data were 
considered. The covariance around the {2, 0.001} point is such 
that it is no longer possible for us to marginally detect or rule out 
an inverted or normal hierarchy directly. However, as explained 
in Sec l3.4l it is still possible to strongly prefer a hierarchy with 
the simple detection of the sum of the neutrino masses if ^2 ™s ~ 
0.1 eV. 



the constraints of both of these parameters from a combina- 
tion of LSS and CMB experiments. 



3.4 Indirectly probing neutrino hierarchies with 
cosmological data at ~ 0.05 eV. 

We have shown that for a fiducial model where ~ 0.25 eV 
and N,^ = 2 there is a good prospect of a very signifi- 
cant detection of (Secs l3.3TT] and l3.3.2p . and a reason- 
able prospect of a detection of Ni, with future SKA data 
(Sec l3.3.51l . The problem is significantly harder if we wish to 
push the bounds to a neutrino mass of ~ 0.05 eV which 
is the limit one would like to attain given that it is the lower 
limit of the parameter space allowed from particle physics 
experiments. 

We have run MCMC chains that assumed 
a fiducial model {f26,f2c,Sl,y,/i,?is,(T8, A^'^} — 
{0.04,0.259,0.001,0.72,1.0,0.9,2}, together with addi- 
tional nuisance parameters, in order to test whether the 
predictions we made hold and to know whether future 
surveys will be able to push the limits for m^. 

We find that for an SKA survey alone the relevant pa- 
rameters will not be measured accurately enough to provide 
strong bounds on the neutrino mass. We plot in Fig |13l the 
MCMC samples from such a theoretical experiment. We find 
that a LSS (SKA) survey on its own will not be able to find 
a signal either from the neutrino density parameter or from 
the number density of massive neutrinos if ~ 0.05. 

It turns out that the parameter degeneracy is too strong 
for there to be a clear detection. A 3a upper limit for the 
neutrino mass would in this case be ~ 0.25 eV (95% 
confidence). 

However the addition of CMB data does improve signif- 



icantly the estimate of at these low masses. As we can 
see from Fig |14l it would be possible to measure the neutrino 
mass accurately with a{m,j) ~ 0.015 if we include both fu- 
ture LSS (SKA) and CMB (Planck) data. This would mean 
that the entire parameter space that is currently still allowed 
by current cosmological and laboratory experiments would 
be probed and a signal from a massive neutrino would be 
detected. 

It is possible for us to ask also whether it is possible to 
detect A'^^ with this signal. We plot in Fig |15l the posterior 
from the MCMC samples for the parameters fi^ and A'^. 
Even though the mass density of neutrino can be measured 
accurately, A^^ is totally unconstrained within the priors we 
have chosen. It would of course with real data make sense 
to then widen the prior and obtain a sensible error for this 
parameter. However, with a strong detection of the neutrino 
mass it is still possible within a certain range to obtain in- 
direct conclusions with regards to the neutrino hierarchy. If 
an inverted hierarchy is assume then we can simply assume 
that there are two massive neutrino and ignore the mass 
of the third neutrino and consider it massless. In this case, 
given that there is a minimum mass splitting, and we find 
a total sum of neutrino masses equal to less than twice this 
value then the hierarchy cannot be inverted. This 'reductio 
absurdum' method implies that if J]] 0.12 eV then the 
hierarchy must be normal. 

We conclude that a combined data set comprising LSS 
(SKA) and CMB (Planck) would be able to detect the signal 
of a massive neutrino any standard scenario. Furthermore, 
if the sum of neutrino masses is below 0.1 eV or above 0.25 
eV there is a good possibility that cosmological data will be 
able to prefer one hierarchy strongly. However if the sum of 
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Figure 13. Illustration of MCMC samples assuming a total energy density in relativistic particles = 0.001 and Nu = 2 which 
corresponds to an absolute neutrino mass scale of 0.05 eV. We consider LSS (SKA) simulated data alone. The vertical red dashed lines 
in the right panel correspond to one, two and three sigma confidence levels. We can see that the parameters A^^ and are not well 
measured, even with exquisite data coming from future galaxy surveys. They alone will not provide evidence for a massive neutrino given 
the lower limit currently allowed by particle physics experiments. 



neutrino masses is in the difficult range 0.1 - 0.25 eV then 
the data we have considered may not be enough for a strong 
preference to be found with regards to the hierarchy. 



4 CONCLUDING REMARKS 

We have shown that future cosmological experiments will be 
able to detect a signal from a cosmological massive neutrino. 
This signal would be encoded as a gentle curvature in the 
galaxy power spectrum correlated with a shift of the CMB 
temperature and polarisation peaks. We argue that it will be 
possible to remove any systematics present in data arising 
from the clustering properties of galaxies by the use of dif- 
ferent galaxy types in different redshift bins. Furthermore, 
data will be available from galaxies with different spatial 
clustering and hence any systematic effects will be assessed 
by comparing data sets with galaxies with different cluster- 
ing and redshift distributions. 

We argue that for a clear detection to be made with 
confidence we need to be sure that any prior used does not 
affect the results significantly. We find that future cosmolog- 
ical data, both LSS (SKA) and CMB (Planck) data will be 
accurate and wealthy enough so that priors will play less of a 
role given that it will be able to measure several key param- 
eters at once including the matter density and the Hubble 
constant. We find as was already pointed out several times 
in the literature, a strong complementarity between CMB 
and LSS data. 

Although it is of utmost importance to measure proper- 
ties of neutrinos directly with particle physics experiments, 
there are several limitations to these techniques given that 
the mass scale of neutrinos is extremely small. Therefore 
other probes from cosmological studies are key to having a 



more complete picture of neutrino physics. Furthermore, if 
we are to assess unusual theories such as scenarios where 
sterile particles are present and have an energy density that 
influences cosmological data, then there needs to be detailed 
comparisons between the values drawn from direct particle 
physics detections and from cosmology. 

We argue here that if the mass scale of the neutrino 
nil, is very small we find ourselves necessarily in a non- 
degenerate scenario. However if this mass scale is large 
enough, the spectrum of masses becomes quasi-degenerate 
and most of the masses are almost equal. Current exper- 
iments cannot yet distinguish between these scenarios but 
may be on the verge of doing so. From a particle physics 
prospective the absolute mass scale of the neutrino is larger 
than 0.05 eV and smaller than 2 eV. 

We have analysed to what extent neutrino parameters 
can be recovered by combining future Planck data and a LSS 
survey, e.g. with the SKA, that would measure redshifts out 
to z ~ 2 over a large fraction of the sky. We conclude that 
this set of experiments would be sensitive over the entire 
mass range allowed by current particle physics experiments 
down to 0.05 eV at the 3(7 level. We find that even though 
the mass of the neutrino can be measured for the entire 
parameter space currently still allowed by experiment, it is 
harder to measure the number density of massive neutrinos. 
We find that down to 0.25 eV it is possible to measure the 
number density well enough to distinguish between a sce- 
nario where one or two of the mass eigenstates are massive 
and hence to be able to distinguish between a normal and 
an inverted hierarchy. 

We have ignored the effects of any varying equation of 
state of dark energy, assuming that this value will mainly be 
determined by the angular scale of the wiggles in the power 
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Figure 14. Wc plot here the equivalent to Fig[7|and FigO using a fiducial cosmology with Qi, = 0.001, and considering both LSS and 
CMB data, alone in the lower plots, and in combination in the upper panel. The vertical red dashed lines in the upper panel correspond 
to one, two and three sigma confidence levels. Wc can see that LSS surveys will, alone, fail to measure a neutrino mass below 0.25 eV, 
and CMB surveys, alone, will also fail to measure such a mass. However, as we can see from the upper pannel, the combination of both 
surveys will be able to measure a neutrino mass as low as 0.05 eV with an error of 0.015 eV. 



spectrum. This information depends strongly on the value 
of w and very weakly on the value of fi^, hence providing 
orthogonal constrain ts as has already been seen in real data 
l|Goobar et al.ll200^ . Of course w = —1 could be the real 
answer. 



We have also chosen to use a single value for Us so 
that the the running of the spectral index rirun is set to 
zero. We argue that the degeneracy that will degrade the 
neutrino sector constraints depends on Us to first order and 
on rirun to second order. The Planck satellite is a mission 
designed to probe the shape of the initial power spectrum to 



unprecedented accuracy and hence will determine Aurun ~ 
0.01 (see Planck 'blue book'Q). • 

We have also assumed the initial conditions to be adia- 
batic pertubations. Although it has been shown that isocur- 
vature perturbati ons may generate degen eracies with neu- 
trino parameters (|Zunckel fc Ferreiralliooil ) . polarisation ex- 
periments from Planck will also impose stringent limits on 
the amount of initial isocurvature. We therefore stress the 
importance of a high quality CMB mission as well as a deep 
all-sky redshift survey in order to determine the neutrino 



^ http:/ /www.rssd.esa.int/Planck 
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parameters, because if this is not the case, limits from LSS 
alone would be also plagued by degeneracies. 

We find that LSS (SKA) and CMB (Planck) data will 
constrain rrii down to 0.05 eV with an error of around 
0.015 eV. If ^ rrii is indeed below 0.1 eV it will be possible 
to rule out an inverted hierarchy indirectly given that having 
an inverted hierarchy with such small ^ mi is in contradic- 
tion with atmospheric neutrino experiments. We therefore 
conclude that if the neutrino mass lies between 0.05 and 
0.1 eV or above 0.25 eV it should be possible to distinguish 
clearly between the normal or inverted hierarchies. However, 
it will not be possible to do so clearly if the mass lies be- 
tween 0.1 and 0.25 eV. These results may be modified and 
improved, if a large range in data points can be used, and 
specially if higher k values can be included in the analysis 
by virtue of improved modelling of non-linear fluctuations. 

It is of great importance that cosmological experiments 
provide us with the same answers and parameters as particle 
physics experiments. Put another way, it may eventually be 
possible to measure any exotic sterile particles that would 
contribute to the cosmic background of relativistic particles 
and hence influence P{k) but which could not be detected 
by particle physics. We argue that it is likely to be possible 
to constrain strongly any of these scenarios given that the 
data will be sensitive down to rrii ~ 0.05 eV and ~ 1. 
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